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Black hole entropy is derived from a sum over boundary 
states. The boundary states are labeled by energy and mo- 
mentum surface densities, and parametrized by the boundary 
metric. The sum over state labels is expressed as a functional 
integral with measure determined by the density of states. 
The sum over metrics is expressed as a functional integral 
with measure determined by the universal expression for the 
inverse temperature gradient at the horizon. The analysis ap- 
plies to any stationary, nonextreme black hole in any theory 
of gravitational and matter fields. 
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Many researchers have suggested that black hole en- 
tropy arises from a sum over boundary (or horizon or 
surface or edge) states . The formal analysis presented 
in this paper shows that this conjecture is essentially cor- 
rect. I consider a system of gravitational and matter 
fields with outer and inner boundary elements, denoted 
B and B' respectively. The entropy for stationary, nonex- 
treme black holes is derived by summing over all states 
on the inner boundary element. The states on B' are 
labeled by the energy surface density e' and momentum 
surface density j'^ pi, and parametrized by the surface 
metric ct^j. The measure for the integration over e' and 
jl is determined by the density of states, which is calcu- 
lated formally as a path integral. Integration over e' and 
J ■ imposes the condition that B' is the bifurcation surface 
of a stationary black hole horizon. The measure for the 
integration over cr^^ is not fixed by the density of states 
and requires further physical input for its determination. 
The needed physical input is the inverse temperature gra- 
dient at the horizon. The inverse temperature gradient is 
a universal quantity, in the sense that it is independent 
of the details of the black hole and independent of the 
theory of gravity under consideration. In this regard, its 
input into the calculation is analogous to the identifica- 
tion of inverse temperature with imaginary time period. 

Recent results indicate that black hole entropy can be 
understood in terms of D-brane states in string theory 
1^ . Although the string theory description of black hole 
entropy is quite compelling, it is limited in scope to a par- 
ticular theory of gravity (as defined by the low energy 
limit of string theory). On the other hand, path inte- 
gral and variational methods have shown that black hole 
entropy appears in all diffeomorphism invariant theories 
of gravitational and matter fields |^,||. This suggests 
that a deeper understanding of black hole entropy can 
be found from within the path integral formalism itself, 
without specialization to a particular theory of gravity. 



The analysis contained in this paper points to such a 
deeper understanding by providing an interpretation of 
the path integral derivation of black hole entropy as a 
sum over boundary states. 

Let me begin by presenting a quick overview of certain 
key aspects of statistical mechanics. Consider a nongrav- 
itating, quantum mechanical system with Hamiltonian 
operator Hy which depends on, say, volume V. (For 
example, the electromagnetic field in a spherical con- 
tainer of volume V.) The partition function is defined 
by Zv{P) — tre^^^^, where (3 is the inverse tempera- 
ture. Eigenstates of Hy are labeled by the corresponding 
eigenvalue, namely, energy E, and are parametrized by 
volume V. The partition function can be written as 



ZviP) = / dEDv{E)e 



(1) 



where Dv{E) is the density of energy eigenstates. That 
is, Dy{E) dE is the number of states with volume V and 
energy in the range dE. An alternative proceedure is to 
define Zv{/3) by a path integral over histories that are 
periodic in imaginary time with period (3. The phase in 
the path integral is given by the action, which is related 
to the classical Hamiltonian Hy in the usual way. The 
density of states Dv{E) becomes a path integral over 
histories with fixed energy E, where E is the value of the 
classical Hamiltonian |^ . 

For a given value of (3, the probability that the 
system will be found in a state with volume V and 
energy in the range dE is VDv{E)dE, where V = 
/ZviP)- The expectation value of energy is de- 
fined by {E) = — 9(lnZy(/3))/9/3, pressure is defined by 
fjp = d{\nZY{f3))/dV, and entropy is defined by 

Sv{f}) = - j dEDv{E)V\nV = \nZv{f3)+f3{E) . 

(2) 

This is the entropy of the system in the canonical ensem- 
ble at temperature 1//9. The microcanonical entropy is 
obtained from the infinite temperature limit, (3 — 0, in 
which V is constant. In this limit Eqs. (1) and (2) yield 
Sv{0) — ln(/ dE Dv{E)), which is the logarithm of the 
number of states with volume V. The total number of 
states N , with no restriction on V, is given by 



N : 



J dEdVfi{V) Dv{E) 



(3) 



Here, ^i{V) is a measure for the V integration that must 
be specified by some physical criterion. This unspecified 
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measure also appears in the definition of the partition 
function for the constant pressure ensemble namely, 
Zp{f3) = J dVii{V)e-l^P^ZvW. In fact, one can view N 
as the infinite temperature limit of Zp{j3). 

The density of states factor in Eq. (1) can be written 
as e'"^^'^*-^-'-' , where for simplicity an arbitrary dimen- 
sionful constant (such as the Planck energy) has been 
omitted from the logarithm. Let us evaluate the inte- 
gral for Zv{(}) in the steepest descents approximation. 
With E = E* {(3) denoting the solution of the extremum 
condition f3 = d{\n DY{E))/dE, we find 

\nZv{l3)^-l3E* +\n[Dv{E*)A] (4) 

where 27rA-2 = _ [d^{\nDviE))/diEy]\j^,. The steep- 
est descents approximation also yields (E) w E* . By dif- 
ferentiating the relation /3 w [9(lnZy(_E))/9i?]|^^^ with 

respect to /3, we find that is given by the mean square 
deviation in energy, A^/(27r) « {{E — {E)Y). A compar- 
ison of Eqs. (2) and (4) shows that 

Sv{P)-HDv{E*)A]. (5) 

Thus, the entropy in the canonical ensemble at tem- 
perature 1//3 is the logarithm of the number of states 
in an energy interval given by the mean square de- 
viation in energy. Equation (5) can be written as 
Sv{P) ~ ln[Zy(£^)A] where [3 and E are related by 
(3^d{\iiDv\E))/dE. 

Now consider vacuum Einstein gravity on a spatial 
manifold S. The boundary dT, need not be simply 
connected — later we will specialize to the case in which 
9S consists of an inner boundary element B' and an outer 
boundary element B. The partition function is defined 
as a path integral on the manifold S x S*^ with actiorj^ 

S"^ =ij dt(^J d^xP'^kj - , (6a) 

H„= I d^x{NH + V'n{) 

+ [ d^x^{Ne - V'j^) , (6b) 

where hij and P'^ are the canonical variables, N and 
are the lapse function and shift vector, and Ti. and Tii are 
the Hamiltonian and momentum constraints. Also, (Tab is 
the metric on 9S and e = fc/(87r) and ji = —2Pijn^/Vh 
are the energy surface density and momentum surface 



'^See Ref. Q. The action in Eq. (6) is defined by the ex- 
ponential in the path integral, and therefore equals i times 
the "Lorentzian action". Also, for simplicity, the condition 
V'Tii = on dT, has been imposed; this means that the 
boundary is at rest with respect to the constant time slices. 
The action in the presence of moving boundaries is given in 
Ref. §. 



density Here, k = (T°'^kab and denote, respec- 
tively, the trace of the extrinsic curvature and the out- 
ward pointing unit normal of (9E embedded in S. 
The variation of the action (6) is 

= i [ dt ( d^x ^{e5N - j,SV' 

J JdT, 

^{Ns-''/2)Saab} , (7) 

where the unwritten terms yield the classical equations of 
motion and (iVs'^V2) = [Nk'''' + {n^d.,N - Nk) a"'']/ {167r) 
defines the spatial stress. The form of SS^-^ shows that 
CTah and the boundary values of N and V are held fixed in 
the variational principle. The partition function, defined 
by the path integral, is a functional of these quantities. 
The boundary value of N determines the period in imag- 
inary time, which is identified with the inverse temper- 
ature on the system boundary: /3 = i J dtN\g-^. (Note 
that /3 is not typically constant on 9S due to gravita- 
tional redshifting and blueshifting.) The boundary value 
of defines the proper velocity of the spatial coor- 
dinate system through the relation f3v'^ — i J dtV^\g^. 
(In the classical approximation, is the spatial velocity 
of the system with respect to observers at rest on the 
boundary |§,||.) The partition function will be denoted 
Zcr[f3,v]. To be precise, Za[f3,v] is the grand canonical 
partition function and plays the role of the "chemical 
potential" for ji. 

Because the Hamiltonian and momentum constraints 
vanish on shell, the value of the Hamiltonian H^r for var- 
ious choices of lapse and shift is given by the densitized 
energy surface density y^e and the densitized momen- 
tum surface density Vo-Ji- These quantities play the role 
of energy. The fact that there are many such "energies" 
refiects the many-fingered character of time in gravita- 
tional physics. The fact that they are surface quantities 
refiects the fact that the gravitational field on 9E is de- 
termined by the mass-energy throughout S. By analogy 
with the nongravitating system discussed above, we are 
lead to the formal identification of ^/ae and y/aji as la- 
bels for (boundary) states. These states are parametrized 
by the boundary metric aab- In direct analogy with 
Eq. (1), the partition function can be written as 

ZAP,v] = J dedjDa[e,j]e-Io^'^''^^^'-^''''^^ , (8) 

where J de dj denotes functional integration over ^/a^^ 
and \foii. In Eq. (8), Dcr[e,j] is the density of boundary 
states; that is, Da[s, j]de dj is the number of boundary 
states with metric aab and energy and momentum surface 
densities in the range dedj. 

The entropy of the gravitational field in the canonical 
ensemble with boundary temperature 1/(3 can be com- 
puted from the partition function Za- [P, v] using the ana- 
log of Eq. (2). Alternatively, one can use the density of 
states I?cr[e, j] and the analog of Eq. (5). In the zero- 
loop approximation, in which only the dominant expo- 
nential contributions to the path integrals for Za[l3,v\ 
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and Dcr[e,j] are kept, the second option is easier. In- 
deed, in this approximation the mean square deviation 
factor A can be dropped and the analog of Eq. (5) be- 
comes Sa[P,v] ~ \n{Da-[e, j]). Here, it is understood that 
(3 — S{lnDc[e,j])/ 6(^/^6) with a similar relation for /3u'. 
Now, from the path integral for Zcr[s,j], we find that the 
density of states is expressed as a path integral over the 
manifold S x S*^ with action 

S'^' = i j dt j Sx{P'^hij -NH- Vn^) . (9) 

In the zero-loop approximation, D equals the ex- 
ponential of the action 5'''^ evaluated at the stationary 
classical solution (if one exists) with the given values of 
e, ji, and dab on 5E. Since S"^'* = for a stationary clas- 
sical solution, the zero-loop contribution to Da[£,i] is 1 
and the entropy Sa[(3,v\ « ln(I?(j[e, j]) vanishes. 

The conclusion Sc\(3,v\ k, is general. Note, however, 
that e and ji label the states for the complete boundary 
dYj of the spatial manifold. If dYj consists of discon- 
nected pieces, then e and ji include labels for states on 
each boundary element. Likewise, (3, u', and tJab spec- 
ify the inverse temperature, velocity, and metric on each 
boundary element. Let us assume that consists of 
an "outer" boundary element B and an "inner" bound- 
ary element B', and make this dependence explicit by 
writing D[e, j;a\£' ^ f \a'] for the density of states and 
Z[l3^v]a\j3' ,v' ]a'] for the partition function. Here, the 
unprimed quantities refer to B and the primed quanti- 
ties refer to B' . The vanishing entropy is expressed as 
S[l3,v]a\(3' ,v'-a'] « 0. 

As an example, consider a Schwarzschild black hole. 
Let the data e, ji, Gab on B and the data e', ][, (t^^ 
on B' coincide with those for two "r = const" surfaces. 
By construction, the classical approximation of the sys- 
tem will consist of the portion of the black hole between 
the two "r = const" surfaces. According to the analy- 
sis above, the entropy of this system vanishes. This is a 
rather intriguing result when one considers the fact that 
the boundary elements B and B' need not lie in the same 
wedge of the Kruskal diagram j|] . If B lies in, say, the 
right wedge and B' lies in the left wedge, the system (in 
the classical approximation) will contain a black hole, yet 
the entropy will vanish. (Note that, in order for the ac- 
tion (9) and the path integral to be well defined, it is not 
necessary for the spatial slices of the classical approxi- 
mation to be nonintersecting.) This, of course, is not the 
usual result for black hole entropy. In the usual path in- 
tegral derivation of black hole entropy only the data on a 
single outer boundary element are specified Loosely 
speaking, the information present in the inner boundary 
data has eliminated the black hole entropy. 

Consider, then, our system with inner and outer 
boundary elements, but with data specified only on the 
outer boundary element. Let us assume that the outer 
boundary data are chosen such that the system is clas- 
sically approximated by a stationary, nonextreme black 



hole. The entropy of the system is the logarithm of the 
number of outer boundary states with metric aab and 
with e and ji in an interval A. The number of such 
states is given by the "density of outer boundary states" , 

D[e,r,n]^ J ds'dj'da' fi[a']D[e,j;a\e',j';a'] , (10) 

obtained by summing over all inner boundary states. 
Equation (8) shows that the integration over state la- 
bels e' and j'^ can be viewed as a transformation to the 
canonical ensemble with /?' = and — 0. (One can 
show that v'^ = as well Q.) Thus, the boundary el- 
ement B' becomes an infinite temperature surface. For 
the black hole that approximates the system, B' is the 
bifurcation surface of the Killing horizon. 

In Eq. (10), which is the analog of Eq. (3) with re- 
spect to the inner boundary data, a measure factor /i[cr'] 
appears. This measure can be specified as follows. Con- 
sider the path integral for D[s,j;a\e',j';a'], constructed 
from the action S"'*. In the steepest descents approxi- 
mation, the integrals in Eq. (10) yield the condition that 
the variation S{S'^^ + ln^[a']) with respect to e' , j^, and 
(T^j should vanish. Equations (6), (7), and (9) show that 
S{S'^'' + ln/i[(T']) includes the inner boundary terms 

^ J dt J {NSiV^e) - V'Si^/^j^) 

+ ^{Ns°-''/2)5(Jab] + I d^x{6\n^i/5aab)S<7ab ■ (11) 

JB' 

From the first two terms we recover the results /?' = 
and = 0. From the remaining terms, and the 

definitions of [Ns""^ /2) and /3, we obtain the condition 
^/a{—n^diP)a°-'' « 167r((51n/i[fT]/(5crab). Now, the factor 
—n^diP, which is the inverse temperature gradient at B' , 
obeys —n^di(3 = 2tt. This is a universal result, in the 
sense that it applies to all bifurcate Killing horizons in all 
theories of gravity. It can be derived, for example, by con- 
sidering the response of a particle detector that maintains 
a constant proper distance from the horizon.^ Derived in 
this way, it is apparent that the result does not depend 



^The boundary data are stationary and nonsingular, so I will 
assume that the appropriate quantum state is the stationary 
Hadamard vacuum ||lO|. Note that such a state exists only if 
the Killing vector field is everywhere timelike. For stationary, 
asymptotically flat black holes such as the Kerr solution, one 
can possibly define the Hadamard vacuum by taking the outer 
boundary B to lie inside the speed-of-light surface. Also 
note that —n^diP equals itiP, where k is the surface gravity 
and P is the coordinate time period [Q. If the black hole is 
asymptotically fiat, and if the lapse function is normalized to 
unity at infinity, then iP equals the inverse temperature /3oo 
at infinity. Under these conditions the result —ri^diP = 2n is 
equivalent to /3oo = 2n/K.. 
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on the particular theory of gravity under consideration. 
Also it does not depend on the details of the spacetime 
geometry since, at sufficiently close distances, any bifur- 
cate Killing horizon is physically indistinguishable from 
a "Rindler horizon" (the horizon in flat spacetime gener- 
ated by a Lorentz boost). To be specific, then, the desired 
result can be derived from the expression P — 27r/a for 
the inverse temperature measured by a detector under- 
going constant acceleration a in the Minkowski vacuum 
of flat spacetime It is not difficult to show that 

1 /a measures proper radial distance from the bifurcation 
surface. Hence, we have —n^di/3 — 2tt and the condition 
on the measure fj.[a] becomes S\n^[a]/6aab ~ ^/aa°-^ 
Integrating this expression, we find that to leading order 
In /i[(T'] « Jg, (Px^/A. 

With the measure /i[cr'] and the previous result 
D[£,j\(j\e',j'\a'] sa 1, the density of outer boundary 
states (10) becomes D[e,j;tT] « exp(A/4). Here, A = 
jg, cPxy/a is the horizon area of the black hole that ap- 
proximates the system. The entropy is then S[(3^v](7] w 
hiD[e,j]a] ~ A/A as expected. 

The above reasoning can be used to derive the black 
hole entropy in any diffeomorphism invariant theory of 
gravitational and matter fields. Start with the general 
action S of Refs. 1^,1). S is the integral over spacetime 
M of a, Lagrangian density C. Assume that a boundary 
term can be added to S to yield an action S^^ whose fixed 
boundary data include the induced metric on dA4. The 
path integral with action S^^ is the partition function. 
As in Eq. (7), the variation of S''^-^ defines the state labels 
e, ji, as well as the spatial stress s°^. As emphasized in 
Ref. , the action S'^^ for the density of states contains 
no boundary terms when written in Hamiltonian form. 
Thus, SS'^'^ contains no variations of or V on d'S, and 
gds jj^^g^ differ from S^^ by boundary terms that change 
the fixed boundary data from A^, y , and aab to e, ji, and 
aab- In particular, SS'^^ must include boundary terms of 
the same form as those that appear in the first integral 
of Eq. (11). 

Black hole entropy is derived by summing over inner 
boundary states to form the density of outer boundary 
states, as in Eq. (10). Integration over e' and implies 
that the inner boundary element B' is the bifurcation 
surface of a black hole horizon. Integration over other 
state labels implies that the products of /3 with certain 
"intensive" variables should vanish at B' Q. The mea- 
sure h[(t' , . . .] depends on the surface metric cr^^ and other 
parameters. It is obtained, as in Eq. (11), from the re- 
quirement that S{S'''^ +ln/i) should vanish for variations 
in cr^j. The spatial stress (A^s°^/2) is needed for this 
calculation, and is obtained from the following consider- 
ations. The actions S and S'^^ differ by boundary terms. 
In Ref. (I it is shown that one of these boundary terms 
has integrand —4:y/aniUj^^'''^djN, and the other bound- 
ary terms are linear in the (undifferentiated) lapse A^ 
or shift V^. Here, U^'^''" is the variational derivative of 
C with respect to the Riemann tensor TZ^^pcr- Also, _L 



designates a component in the direction, where is 
orthogonal to E. Now, 5S contains only boundary terms 
at dM. If the lapse on an element of 91], say B' , van- 
ishes, then B' is a surface where the foliation degenerates. 
In this case the only possible contribution at B' to the 
boundary terms of 5S is a "corner" term proportional to 
5{^^N) 0. On the other hand, we see that 5{S - S"*") 

includes a term jg,{~2y'aniU^^''^djN)a"'''6(Tab- This is 
not a corner term, so it must come from SS'''^. There- 
fore the spatial stress, which appears in the bound- 
ary terms of SS'^'^, must be given by {Ns°'^/2) — 
{2niU^^'' dj N)(t'^'' + • • •. Here, the unwritten terms 
contain undifferentiated linear factors of A or y. The 
condition on the measure that follows from Eq. (11) is 

then 2^/aniUQ^-'^nj{—n^dkl3)(7°'^ — 5\ti^[(t,. . .]/6aab, 
where diN is proportional to With the uni- 

versal result —n^dk(3 = 2tt for the inverse tempera- 
ture gradient, we find that to leading order the mea- 
sure is given by ln/x[CT', . • ■1~ Stt /g, ^/aniUQ^''^nj. This 
yields the desired result S[I3,v]cf] « \iiD[e,j]a] ~ 

—2tt /g, -Jae^^Uf^"'^'^ epa, where e^^ — 2u[pn^] is the bi- 
normal of the bifurcation surface B' . 

I would like to thank J.W. York for valuable insights. 
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